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ABSTRACT 


An analytical approach to determine an optimum laminate for a 
variety of thermal and mechanical loading combinations is presented. 
The analysis is performed for a linear elastic material under static 
mechanical and uniform thermal loadings. 

The problem is restricted to a unit width and length laminate 
with angle orientations resulting in an orthotropic, symmetric and 
balanced configuration. This allows for the elimination of bending- 
extention coupling effects, furthermore only in-plane loads and uni- 
form temperature change through the thickness are considered. Thus 
the problem is reduced to that of a plane stress state. 

The optimization is performed by a general purpose program, 

AESOP (Automated Engineering and Scientific Optimization Program) 
developed for NASA by the Boeing Co. An objective function defining 
total strain energy, is formulated and an optimum laminate design 
determined subject to constraints on stiffness, average coefficient of 
thermal expansion, and strength. The objective function is formulated 
in terms of the orientation angles, number of plies and material 
properties. 

The method presented has, in varying degrees, shown that the de- 
sign of a laminate can be accomplished using strain energy minimiza- 
tion as the primary criteria. It is felt that by minimizing strain 
energy, reserve strength is maximized. The inclusion of a failure 
criteria may result in a non-feasible solution based on the purpose of 
a failure criteria to maximize stress, within failure bounds. The 
results of various combinations of applied constraints in the opti- 
mized design process are presented and discussed. 
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Chapter 1 
INTRODUCTION 

1.1 GENERAL 

The goal of structural engineering is* for the most part, to 
arrive at a design which Is the most economical and at the same time 
satisfies all load and deflection requirements. This goal is best met 
by implementing some optimization technique to aid in the design pro- 
cess. The use of optimization is a relatively recent development in 
the design of structures and has not yet enjoyed the intensive study 
needed to yield the potentially powerful tool it promises to be. 

Design, at present, usually amounts to a series of analysis procedures, 
or iterations, until a suitable solution is found, however with the 
advent of modern high speed digital computers the use of an optimiza- 
tion technique to solve the same problem can yield a superior design, 
and a savings in computational effort. 

Laminated composite materials, being composed of many plies at 
various angles of orientation, readily allow the direct application of 
optimization techniques to arrive at an optimum design. Previous 
Investigations have dealt with this problem on the basis cf a minimum 
weight design with various constraints applied. Many of these efforts 
have been of the "try them all" approach. 

The "try them all" approach amounts to the systematic perturbing 
of angular orientations by a set increment until all possible configura 
tions have been analyzed. This is fine for a very small number of 
plies in the laminate but the total computational time involved for 
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a large number of plies is prohibitive. This leads to the use of non 
linear programming techniques to arrive at an optimum design. 

Little attention has been paid to the effects of thermal loading 
on the design, or the need for tailoring laminates for certain thermal 
expansion characteristics. This type of problem is becoming critical 
as higher temperatures environments and extreme temperature ranges are 
encountered, especially in spacecraft and hypersonic flight. 

1 .2 OBJECTIVES 

The primary objective of this report is to develop a computer code 
which can be used to find an optimum composite laminate design, given 

the mechanical and temperature loadings. In addition, the capability to 
use the code purely for analysis is provided. 

The design is found by perturbing the design variables until the 
objective function, consisting of total strain energy, is minimized 
subject to constraints on stiffness, coefficients of thermal expansion 
and failure criteria. The use of the total strain energy of the 
laminate as the objective function to be minimized is thought to be 
a means of maximizing the total reserve strength the laminate has 
available. 



Chapter 2 
LITERATURE REVIEW 

A number of schemes for optimum laminate design have been formu- 
lated during the early years of composite material technology develop- 
ment. The earliest of these deals with fiberglass laminate design by 
Hackman and Stotler [1]. This procedure uses the number of plies, 
individual ply material type and ply orientation as the design 
variables. A combined load envelope is developed and adjusted to 
reflect fatigue and environmental factors. The design is then found, 
by means of a polar loading diagram, through a series of material and 
ply orientation selections. This type of design method lacks the 
ability to consider ply interaction. 

The computer program RC7 written by General Dynamics [2] de- 
termines an optimum laminate by means of a series of ply additions 
and reorientations. RC7 has the capability to handle 20 plies, con- 
tains two strength criteria, Tsai -Hill and maximum strain but stiffness 
is not considered. The ply angles are allowed to vary from 0 to 180 
degrees in increments of 5 degrees. This procedure would seem to be 
very good for their laminates but the large number of reorientations 
involved tend to negate the value for a large number of plies. 

Approaching the problem by a different method, Foye and Baker [3] 
use a random search technique as opposed to the exhaustive, systematic 
search. Initially this procedure determines the number of plies by a 
random search method then in the second step uses non-linear pro- 
gramming techniques to find the optimum orientations. 
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An analytical approach is taken by Bush [4] to the optimization 
problem. This procedure uses a unique simultaneous equation solution 
approach formulated on the basis of the load-stress and stress-strain 
relations. The design variables include total thickness and ± angles. 
The design is constrained to be orthotropic about the principle load 
directions. 

Another type of laminate design optimization procedure has employed 
families of plies. Two such programs, OPLAM by Grumman Aerospace [5] 
and OPTLAM by Douglas Aircraft [6], follow a procedure similar to RC7 
using [0/±45/90] orthotropic laminate. Basically, these two procedures 
go through a series of ply reorientations, restricted to the previously 
mentioned angles, and ply deletions until a minimum ply solution is 
obtained. 

The method found in the Advanced Composites Design Guide [7] is 
a more general constrained minimization program. This procedure uses 
the method of centered circles to arrive at a constrained optimum 
solution. 

Schmit and Farshi [8] present a method for optimum laminate 
design using a series of linear programs to approximate the nonlinear 
programming formulation. The solution implements the method of in- 
scribed hyperspheres to obtain the optimum laminate. This procedure 
considers multiple in-plane loads, strength and stiffness. 

All of the methods reviewed above involve laminate design. The 
next logical step is the assemblage of these laminates into structural 
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elements. In references [9-28] these types of applications are discus- 
sed in great detail. The general structural element application usually 
involves finite element analysis in conjunction with a laminate analysis 
program. A further step in this direction involves the design of com- 
plex structural elements, such as complete wings. The details of these 
applications are discussed in references [29-31] which are beyond the 
scope of this investigation. 

The methods presented in the detailed review, include, systematic 
search techniques, random search techniques, analytical approachs and 
more sophisticated non-linear programming methods. Many of these 
methods tend to be computationally exhaustive for a large number of 
plies, and as a result this investigation shall deal with the usage of 
a series of non linear programming techniques as a means of obtaining 
an optimum solution. 

All of the methods reviewed used weight as the function to be 
minimized with various combinations of constraints, either strength or 
stiffness generally employed but coefficients of thermal expansion 
were neglected for the most part. 

The current investigation shall be restricted to laminate designs 
with thermal effects included. Total strain energy shall be employed 
as the function to be minimized, subject to stiffness, strength, and 
coefficient of thermal expansion constraints. 



Chapter 3 

THEORETICAL FORMULATION 

3.1 LAMINATE ANALYSIS 

The laminate analysis is based on a single, integral structural 
element made up of two or more laminae bonded together/ The indivi- 
dual lamina properties govern the response of the laminate. Logically 
the lamina is the initial step in formulating the laminate analysis. 

The lamina is assumed to be a homogenous orthotropic material in 
a plane stress state. Orthotropic materials have three planes of 
material property symmetry (Fig. la). The general constitutive equa- 
tions for a lamina in its natural coordinate system (Fig. lb) may be 
expressed as: 



11 

S 12 

S 16 

( a l j 


12 

S 22 

S 2 6 

J °2 

(3.1) 

16 

S 26 

S 66- 

U 12 ) 




6 



2 y 



(a) ORTHOGONAL LAMINA (b) NATURAL LAMINA (c/ TRANSFORMED LAMINA 


COORDINATE SYSTEM COORDINATE SYSTEM “ 

FIGURE I. LAMINA COORDINATE SYSTEMS 



8 


The strain-stress relations of (3.1) can by matrix inversion, yield 
the stress-strain relations: 



Q 16 = f? 26 = 0 


The stresses and strains must now be expressed in terms of an 
arbitrary coordinate system (Fig* Tc) rotated by an angle theta from 
the x-axis to the 1-axis. These stresses and strains can be expressed 
in the following form. 



(3.5) 
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and 



where the coordinate transformation matrix, T, is 



~c 2 s 2 2 sc 

s 2 c 2 -2sc 

2 2 

L-sc sc c-s_ 


and the superscript -1 denotes the matrix inverse. 



~c 2 s 2 -2sc " 

2 2 

c c 2 sc 

2 2 

_ sc -SC C 


(3.6) 


(3,7) 


(3.8) 


where 

c> cosine e 
s = sine e 

In (3.6) it is noted that strains transform via the same transforma- 
tion if the tensor definition of shear strain is used (engineering 
shear strain divided by two). However, if the matrix 



' 1 0 O' 

0 1 0 
.0 0 2 . 


(3.9) 
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due to Reuter [32] i$ implemented, the strain vectors become, 


- CR] 


(3.10) 


(3.11) 


Thus, in order to obtain the stress-strain relations in the x,y coordin- 
ate system in terms of the material properties and the laminae orienta* 
tions equations (3.3), (3.5), (3.10), [3.6) and (3.11) are combined to 
yield. 


= [T]” { a 2 ? = [T3‘ , [Q][R][T][R]’ , {e y 


( 3 . 12 ) 


however [R][T][R3 _1 - [T]“ T where superscript T denotes matrix trans- 
pose, thus abbreviating 


[Q] = [t]' 1 [q][t]- t 


(3.13) 


and in simplified form (3.12) becomes. 



n 



The derivation for the strains follows a similar path and can be found 
in more detail in [33] Jones, 

The previously expressed relationships are valid for any lamina 

th 

of a laminate. Thus the stresses for the k ply are. 


{<•• } - [Sllcl (3.15) 

X k K X k 

In order to derive the stress strain relationships for the lamin- 
ate the Kirchhoff hypothesis for plates must be introduced to 
account for the variation of strain through the thickness. This 

±. L. 

representation of strain for the \C U ply in terms of the mid-plane 
strain and curvature may be written as 



Where z is the distance from the laminate midplane to the lamina mid- 
plane and the stresses in the k** 1 ply as. 
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The strains vary linearly through the thickness but due to possible 
variation of the Q matrix from ply to ply the stress variation will, 
in general, be nonlinear. 

The resultant laminate forces (force per unit length) and moments 
(moment per unit length) (Fig. 2) are obtained by integrating the 

j . : ■ 

stresses of each ply through the thickness (Fig. 3) as follows. 


f H n f z k 

(N) =/ {0} dz « E / {a } dz 

J- H x k x k 


* i, 


{a..} zdz 


■\f h 

>-<K, 


{0 } zdz 

A 


(3.18) 


Substituting (3.17) in (3.18) yields. 


n 


{N} = e [QL< 
k=l K 


j j 

Jz k-1 J 


2 k 


(xJzdz 


{M} = e CQ] k |l {e°}zdz + j {K}z c dz 

r z k-l 7 

I K 1 z k-l 


z k-l 


z> 

/ 


(3.19) 


Simplifying (3.19) yields. 


IN} = [A]te°> + [B]{k> 
(M) * [B]{e 0 } + [D]{<} 


(3.20) 


A ij = ^ C 5 ij> k Cz k- Z k-l) 


where 
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FIGURE 2 FORCES AND MOMENTS ON 

A FLAT LAMINATE 
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B ir! j, ^j) k t* 2 k -*k-i) 


(3.21) 


°ij = i ( ^j> k tz k- 2 k-i> for i > i - 1 ' 2 > 6 


Rewriting (3.20) in combined form results in 


(■"*) = [W] K“) 


or in inverted form 


(tI = [b : H^]I'"'I 

Complete details of this derivation may be found in [33] 


(3.22) 


3.2 THERMAL ANALYSIS 

The thermal analysis for laminated composites is developed by 
following the strain history of the individual lamina from the un- 
cured state through curing to the final bonded state. The material 
properties of the lamina are assumed to be constant through the 
temperature range, thus neglecting any degradation effects. Com- 
bining these thermal curing effects with the effects of increased or 
decreased operating temperature on the laminate, the total thermal 
analysis can be developed. Initially, examining the coefficients of 
thermal expansion in the natural coordinate system, the strains due 
to a temperature change for the k ply can be obtained. Thus the 
expressions for the free thermal strains in natural coordinates are. 
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e.j = a-[ AT 


g 2 S a 2 AT 
1.L 

or more simply in vector form for the k u ply are, 



(3.23) 


(3.24) 


The strains of (3.24) may be transformed to x,y coordinate system by 
means of (3.6) and (3.11) as expressed below, 

“1 ) 

“2 / aT= 

0 'k 

where the coefficients of thermal expansion are transformed to the x,y 

coordinate system and noting the a term as an apparent coefficient of 

xy 

thermal shear. The expressions for stress resulting from thermal strain, 
if the lamina is restrained, may be found by substituting (3.24) in (3.3), 
yielding, 



(3.25) 






(3.26) 


Also in x, y coordinates the stresses may be found by substituting 
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(3.25) in (3.14) yielding. 



(3.27) 


The next step is to develop an expression for the coefficients of 
thermal expansion, for the laminate. Substituting (3.27) in (3.18) to 
obtain the equivalent forces due to the thermal loadings for a 


symmetric laminate [B] = 0 yields, 


(N) = AT/ [Q] k (a x } dz 


and 


t N } = [A]{e°> - AT " [Q] k t“ x } k ( z k - z k .i> (3.29) 

k-1 

where {N> is the equivalent thermal force and {e 0 } is the equivalent 
thermal midplane strain which is uniform through the laminate. Thus, 
if (3.29) is rearranged to yield midplane thermal strain, 

<e 0 } • [A]" 1 " CEJ] k {<x x > U k -Z K .,)AT (3.30) 

k- lk 

where, simplifying to separate the laminate coefficient of thermal 
expansion yields, 

n 

{o } = Z [Q] k {« x l ( Z k - Z k -T) (3.31 ) 

k" 1 k 

The analysis via the strain history is now developed. The indi- 
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vidual plies are laid together in the specified angular orientations 
at ambient temperature and are then subjected to pressure and heat. 

The temperature is raised to a maximum and then decreased through the 
stress free temperature down to ambient, the specifics of the cycle 
may vary but this is a general description. The lamina are considered 
cured into a laminate after the stress free temperature is reached and 
no longer act individually. The laminate behaves as a single struc- 
tural component from the stress free temperature down to the ambient 
temperature and then either up or down to the operating temperature. 

A typical laminate cure cycle may be seen in (Fig. 4). Due to this 
change from the unbonded to the bonded state at different temperatures, 
each effect must be accounted for in the analysis. With this tempera- 
ture history in mind the strains of the laminae (Fig. 5) are formulated 


{e v } = {c.} AT 1 + {a} AT 2 + {5} AT 3 
X 1 X, *1 

(3.32) 

<c v > ~ {«„} AT 1 + {«} AT 2 + {a} AT 3 
X 2 * 2 * 2 

and, 

AT 1 + AT 2 + AT 3 =0 (3.33) 

thus, 

AT 3 = - (AT 1 + at 2 ) (3.34) 


The strains produced in the curing process will be designated 
residual thermal strains and expressed as. '. . . 



AT = AT 3 -i- AT* 
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{e y ) R = {».„> (-AT 3 ) + {a} AT 3 = AT 3 ({a} - {«„} ) 

X T X 1 1 

(3.35) 

{e } R = {o y } (-AT 3 ) + (5) AT 3 = AT 3 ({a} - {a y l ) 

X 2 X 2 x 2 

or in general form, 

{e y } R = AT 3 [{5} - {«} ] (3.36) 

x k x k 

where, 

{e°) R = AT 3 {5} (3.37) 

and the equivalent forces 

{N} R = [A]{e°} R (3.38) 

Substituting (3.36) in (3.14) results in the expression for the 
equivalent stresses, 

{0 R = CQ]|.AT 3 ({5} - {a } ) (3.39) 

x k K x k 

In a similar manner the expressions for the stresses and strains due to 
the temperature change from ambient to operating can be derived. 

During this change the laminate is bonded and no effects for the in- 
dividual lamina are included, thus, 

{e„} E = {e°} E = AT 4 {5> (3.40) 

k 

and the additional equivalent forces are, 

{N) E = [A]{e°} E (3.41) 

Substituting (3.40) in (3.14) results in the expression for the 
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additional equivalent stresses 

<»/ “ [Q] k AT 4 (5> (3.42) 

k 

The total temperature effect may now be formulated by combining the 
two effects, however, first the term aT may be defined as, 

AT = AT 3 + AT 4 (3.43) 

The results of (3.43) may be used to simplify the equations that 
follow, thus the total midplane strain is, 

{e°}T = {e°} R + {e°} E = aT{3) (3.44) 

and total equivalent forces are, 

{N} T = {N} R + {N> E = [A]{e°} T (3.45) 

The lamina strains may be expressed as, 

{e ) T = {e } R + {e } E = aT<«> - AT 3 {a> (3.46) 
x k x k x k x k 

and the resultant equivalent stresses are, 

(a y >^ = (a I** + {a } E = [QU . aT ( { ot} - (a } ) (3.47) 

x k . k x k K x k 

Thus, if the cured state is restored (aT= 0) or in other words the 

lamina strains are present but the stress-free state is satisfied as is 

now shown, 

{e v } T = AT 3 {a } (3.48) 

x k x k 
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fo y } T = 0 (3.49) 

X k 


3.3 LAMINATE ENGINEERING CONSTANTS 

The laminate engineering constants may be determined for a 
symmetric laminate by applying a unit inplane load, only in the 
direction of the desired constants. This may be shown by the fol- 
lowing relationship. 



where the average stresses {J } are, 

A 



and 2H is the total thickness. 
Defining Youngs modulus as. 


(3.48) 


(3.49) 


E _!i i=x,y,xy (3.50) 

1 ' e i 

Thus for N x equaling unity the relationships for average stress and 
midplane strain become 



* N„4t) 


(3.51) 


OBI6INAL PAGE ® pooa 
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and 


4 = a-'n, 

x n X 

Substituting (.3.51) and (3.52) in (3.50) yields, 

, . "x^. 'iff 1 

E x " T^Tir " 7T " TH — 


A li N x A ,1 A i i ( 2H ) 


In a similar manner E and 6 are 

y xy 


E.= 


1 


and 


y A^(2H) 


. . 1 
’xy ■ a -i 


Ag 6 ( 2H ) 


(3.52) 


(3.53) 


(3.54) 


(3.55) 


The expressions for Poisson's ratio may be derived by the same 
method, first defining Poisson's ratio as 


v ij = - if 


yields 


xy 


e a’ 1 
= . = . _2! 


and 


Sr 

Si 


yx 


(3.56) 


(3.57) 


(3.58) 


y "22 

Further constants of interest are those of [34] by 
Lekhnitski called coefficients of mutual influence. These are simply 
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presented here* more detail may be found In the previously mentioned 
reference. 



3.4 FAILURE ANALYSIS 

The failure analysis is formulated on the basis of the Tsai -Hill 
[35] criteria for a lamina. Failure of the laminate is based on 
the first ply failure theory such that if one ply fails the entire 
laminate is considered to have failed. This criteria may be expressed 
in terms of the lamina natural coordinate stresses and the individual 
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lamina strengths as follows. 


°1 

~k 

r 


°\°z 


2 

x 12 


.2 - 


< 1 


(3.65) 


This criteria was chosen over the Tsai-Wu theory [36] in order to 
avoid the laboratory determination of a key value needed for the 


evaluation. 

The Tsai-Wu theory is of a more general nature and offers a more 
responsive strength determination due to the greater number of terms. 
This greater accuracy requires the knowledge of an experimental value 
that must be determined by means of a biaxial test to failure. Also, 
the results from the two theories differ only to a small extent over 
the range of most angular orientations. 


3.5 STRAIN ENERGY 

The total strain energy of the laminate is determined by means of 
a summation of the components strain energy due to the various mechan- 
ical and thermal loadings. 

The general expression for the elastic strain energy of a linear, 
elastic body may be found by considering an infinitesimal element 
(Fig. 6 ) of dimensions dx, dy and dz. Thus for the uniaxial stress 
state the force acting on the right or left face is a dy dz where dy 
dz represent an infinitesimal area of the element. Due to this force 

the element deforms an amount e dx where e is the strain in the x 

x x 

direction. Having assumed linear, elastic material, stress is pro- 
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rIGURE 6. ELEMENT IN UNIAXIAL TENSION 
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portions! to strain and further, the force acting on the element 
Increases from zero to its full value In a linear value. The average 
force acting on the element during the time that deformation is taking 
place is a dy dz/ ? . The average force multiplied by the distance 
through which it acts is the total work done on the element. There 
is no dissipation of energy for a perfectly elastic body and therefore 
the work done on the element Is stored as recoverable internal strain 
energy. The strain energy may be expressed as 


dU-- (id dy dz) • (e dx) = i a e v d V (3.66) 

d x X C X X 

where dV is the element volume. Rearranging (3.66) yields the strain- 
energy density 


.. | 0 £ 

dll _ n _ X X 

dV " u " 2 


(3.67) 


Similar expressions may be derived for Oy and x X y. 

The total strain energy for the laminate is expressed as. 


i n MM T T 

u = l Z ({a.} {£.} + {a.} {e.} ) (3,68) 

c k=l 1 k 1 k 1 k 1 k 

where repeated subscripts i indicate summation and M denotes mechani- 
cal loading effects. 



Chapter 4 

OPTIMIZATION FORMULATION 

4.1 GENERAL 

The optimization procedure used Is the AESOP program described in 
Ref. [37,38]. Basically, this program is a series of multivariate search 
techniques for non-linear systems. The optimizer is easily coupled to 
the synthesis program by means of storage linkages. The synthesis 
program computes the objective function and the constraints functions 
which are then supplied to the optimizer for evaluation. The optimizer 
then perturbs the design variables until an optimum design, consistant 
with constraint conditions, is obtained. Nine different search tech- 
niques may be employed, either separately or in any combination to seek 
an optimum, thus allowing freedom from method dependent solution pro- 
blems. A maximum of one hundred design variables are permitted while 
as many as twenty constraints functions may be utilized. Details may 
be obtained from the references. 

The search techniques available in AESOP and a brief description of 
each follow, 

(1) Sectioning - series of one-dimensional searches parallel 
to the coordinate axes. 

(2) Pattern - search in direction of the previous favorable 
search. 

(3) Magnification - search in direction of gain due to pro- 
portional change at all parameters . 

(4) Steepest-Descent - search along the weighted gradient 
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direction 

(5) Adaptive Creeping - similar to (1), but small increments in 
general minimum direction. 

(6) Quadratic - similar to (4) but search direction is along a 
sequence of second-order surfaces as opposed to first-order 
surfaces of (4). 

(7) Davidon's Method - series of searches of (4) type to 
approximate a search of type (6). 

(8) Random Point - evaluation of function at a set of uniformly 
random points. 

(9) Random Ray - search along a sequence of random rays having 
uniform distribution. 

4.2 MATHEMATICAL FORMULATION 

The mathematical formulation of the multivariate optimization is 
presented for the general case without any technique dependent aspects. 
Basically, the goal is minimization of an objective function of the form, 

$ = $(B.) i=1 ,2, . . . ,m (4.1) 

where 3^. are the design variables subject to a system of constraints, 

Tj = j=i ,2, . . . ,P (4.2) 

Constraints may also be applied directly to the independent variables 
by specifying a feasible control space as follows, 

6j 5 8, < e" 1=1,2 m (4.3) 
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Further, constraint functions of an equality nature may be treated 
as an unconstrained problem by replacing the objective function with a 
penalized objective function as follows 

P o 

$ = 0 + z W. * ¥4 (4.4) 

j=i J 3 

where if the W. (the positive error weighting constants) are sufficient- 

J 

ly large in magnitude, minimization of (4.4) is the equivalent of mini- 
mization of (4.1) subject to (4.2). 

This method is referred to as the exterior penalty function ap- 
proach. Use of this approach results in the design only being valid 
if all constraints are satisfied, thus only the final solution, if met, 
is a valid design. This is opposed to an interior penalty function 
approach where each design satisfies a set inequality constraints and 
each successive design approaches the minimum. 



Chapter 5 

RESULTS AND DISCUSSION 

The computer program developed in accordance with Chapter 3 and 4 
was used to solve a series of example problems to demonstrate the 
various design options. These example problems used Boron/Epoxy, AVCO 
5505/4 as the design material. 

Designs were determined for one mechanical and one temperature 
loading case with constraints on stiffness and/or coefficients of 
thermal expansion and/or ply strength. These desired constraint 
values correspond to the material properties of titanium, aluminum and 
an all composite [0,±45,0] 8 ply, symmetric, Boron/Epoxy laminate 
(Table 1). 

The mechanical loads represent a case where ply failure will be 
critical and therefore tax the strength criteria. Temperature loads 
were chosen for a typical cure and operating environment. These 
loads will be used for all problems unless noted otherwise (Table 2). 

Additionally, all examples were run for 300 function evaluations 
with all other optimization parameters being identical. 

The design options are broken down into the various catagories in 
the results. The major division is whether or not the strength cri- 
teria is considered as in sections 5.1 and 5.2. These sections each 
have subsections according to which method is chosen, either letting 
the angles vary from -90° to 90° for a set number of plies or to let 
the number of plies vary for each set ply angle. Only the first of 
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TABLE 1 

MATERIAL PROPERTIES 



Boron/Epoxy 

Boron/Epoxy 

Titanium 

Aluminum 


AVCO 5505/4 

CO, ±45,0] 



E v (msi ) 

30.0 

17.05 






15.8 

10.4 

E y (msi ) 

2.7 

5.23 



6 xy (msi) 

0.93 

4.43 

6.0 

3.9 

v 

0.21 

0.67 

0.34 

0.33 

xy 





a y (v1n/in/°F) 

2.5 

2.28 



X 


l 

4.8 

13.1 

a y (viin/in/°F) 

13.1 

6.45 



Xj(ksi } 

188 




X r (ksi ) 

-362 







140 

80 

Y T (ksi) 

9.1 




Y c (ksi) 

-45 




Z(ksi ) 

19.2 
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TABLE 2 

MECHANICAL AND TEMPERATURE LOADS 


N = 1440 Tb/in 

A 

N = 1440 Tb/in 

H = 1020 Tb/in 

xy 

AT 3 = -200° F 
AT 4 = 100° F 
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these options is available in section 5.1 for design without a strength 
criteria. It is felt that the problem of constraints on material 
properties, or material property matching, can best be resolved by means 
of reorientation alone. Furthermore, the option for perturbing the 
number of plies is felt to be needed primarily for strength considera- 
tions. It may be noted that for the option of perturbing ply angles 
for a set number of plies that the number of plies will fluctuate. This 
is due to the fact that for angles of the range -90° < e < 0° and 
0° < 0 < 90° (angle ply laminate) a positive theta ply is always ac- 
companied by a negative ply . 

Initially, a test problem to check out the minimization process was 
run for temperature loading only and without constraints. Intuitively, 
it was felt that this should yield some unidirectional laminate, having 
a total strain energy of zero. Three problems, using 2, 4 and 8 plies, 
yielded the following results (Table 3). The results demonstrated that 
the computer program was performing as expected and on this basis the 
following examples were run. 

These designs are constrained to match various, feasible material 
properties while minimizing the total laminate strain energy. The 
results obtained are presented and discussed to validate the design 
process employed. 

5.1 DESIGNS WITH MATERIAL PROPERTIES CONSTRAINED 

This section deals with the problem of using a composite laminate 
to match the material properties of another material. There are many 



TABLE 3 

UNCONSTRAINED STRAIN ENERGY RESULTS 


Number of 
Plies 

Final 

Configuration 

Strain 

Energy 

Number of 
Evaluations 

2 

[24] s 

0.3980 E-14 

17 

4 

[7.4,7.6,8.0,7.2] 

0.9048 E-03 

68 

8 

[R.8,7.6,8.0,8.0] ; 

0.3873 E-02 

134 
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occasions where strength Is not of primary consideration, due to the 
fact that the design Is minimum thickness critical, such as aircraft 
skin designs. Another area where material property considerations are 
critical is that of extreme, thermal deflection requirements. 

These applications represent a large part of the present thrust 
of composite hardware technology, which deals with the replacement of 
present hardware with a composite design. The importance of property 
matching is illustrated when the composite replacement part is compared 
to the present part, for properties being equal, the composite is 
advantageous on the basis of weight. 

5.1.1 DESIGNED BY PERTURBING PLY ANGLES FOR A SET NUMBER OF PLIES 

This approach to a non-strength critical design relies entirely 
on the angles or orientation to satisfy the stated constraints and 
locate a minimum strain energy laminate. The possible number of com- 
binations of angles increases as the number of plies increase, there- 
fore the number of plies chosen for the starting point deserves serious 
consideration and should be adequate to allow the process the proper 
amount of freedom to satisfy the various constraints. The examples of 
this section were run for 2 plies, thus the maximum number of plies is 
8 for an angle ply laminate and 4 for a unidirectional or cross ply 
laminate. The factor of two is due to the assumption of symmetry for 
the laminate. ' - 
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5. 1.1.1 COMPOSITE LAMINATE TAILORED TO MATCH THE PROPERTIES OF A 
[0,±45,0] LAMINATE 

The use of an all composite example is to check out the program for 
convergence to a known feasible design. The boron/epoxy [0,±45,0] was 

i 5 

analyzed, using the analysis option, to obtain the laminate properties 
that were given earlier in this chapter. Seven various cases were run 
for this target material for the standard mechanical and temperature 
loads, unless noted. The constraints and nonstandard conditions for 
these were, as follows; 

(1) E , E and 6 constrained; 

x ) xy 

(2) E and E constrained; 

x y 

(3) a and a constrained; 

x y 

(4) E , Z and G constrained, N = 0; 

a y ajt y 

(5) E x , E y and G xy constrained, N y = N xy = 0; 

(6) E x , E y and G xy constrained, N x = N y = N xy = 0; 

(7) E y , E and G constrained, aT 3 = Ar = 0. 

x y xy 

These results are shown in Table 4. 

Cases (1), (2) and (3) offer the best overall convergence to target 
properties and also demonstrate reduced strain energy at between 5 to 
10 percent compared to the target material. Case (3), oddly enough, has 
good convergence to properties other than those constrained and has 
excellent thermal tailoring ability. 

The nonstandard loading conditions yield reasonable property 
tailoring, with (4) being the most responsive for variations of 
mechanical loads. Deletion of temperature loads, in example (7), 
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TABLE 4 


COMPOSITE LAMINATE - {0,145,03 RE5ULTS FOR PLY ANGLES PERTURBED KITH MATERIAL PROPERTIES CONSTRAINED 




1 

5 

2 

3 

4 

Cases 

5 

6 

7 


Target 

•v E y a " d G xy 

E x and E y 

a and a E 

x y 

' , E and G 
x y 

E , E and G 
xy V y xy 

E , E and G 
x y xy 

E , E„ and G, 
x r 


Material 

Constained 

Constrained 

Constrained 

Constrained 

v° 

Constrained 

vv° 

Constrained 

N =N -N =0 
x y xy 

Constrained 

aT 3 =4T 4 =0 

Laminate 

[0,l45,0] s 

[i5I,±7] s 

[i50,±7] s 

[±45, t 12] s 

[*50,t4] s 

rf20,i25] s 

[+14.i30] s 

[s« ) +S] s 

E (msl ) 

17.05 

15.87 

15.92 

15.05 

16.35 

•’7.51 

18.07 

16.30 

x t 


-6.9 

-6.6 

-11.7 

-4.1 

+2.7 

+6.0 

-4.4 

E {nisi 5 

5.23 

6.70 

6.41 

4.92 

6.49 

2.53 

2.85 

5.19 

y X 


+28.2 

+22.6 

-5.9 

♦24.0 

-51.6 

-45.5 

+18.4 

G (msi) 

4.43 

4.48 

4.53 

5.01 

4.39 

4.43 

4.32 

4.46 

*y j 


+1 .2 

+2.3 

+13.0 

-0.8 

0.0 

-2.5 

♦0.7 

V 

0.67 

0.54 

0.57 

0.76 

0.55 

1.25 

1.12 

0.58 

xy 

t 


-19.4 

-10.9 

+13.4 

-17.9 

+86.6 

+67.2 

-13.4 

a /-in/in/°F) 

2.28 

2.59 

2.53 

2.17 

2.54 

1.23 

1.47 

2.43 

V 


+13.6 

+11,0 

-5.0 

+11.4 

-46.1 

-35.5 

+8.8 

-a (u'n/in/ 0 F) 

6.45 

5.30 

5.46 

6.47 

5.47 

11.85 

10.80 

5.64 

y 


-17.8 

-15.3 

+0.3 

-15.2 

+83.7 

+67.4 

-12.6 

Strain energy 
loading) 

(varies with 

1356 

1369 

1424 

841 

297 

9 

1267 

.* 

NOTE: 

-9.7 

all percentages 

-8.9 -5.2 

refer to variations with respect 

-3.7 -10.0 

to the target iraterial. 

-51.4 

-7.1 
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also resulted in rather good overall property convergence. Cases (5) 

and (6) exhibit a fairly wide spread of values other than E and G , 

x a y 

The larger variations of E can be attributed to its smaller value, 

♦/ 

relative to E and G combined with use of equal error constraint 

A A Jr 

weighting factors. 

Using nonstandard load conditions demonstrates the ability of the 

optimization procedure to respond to different load environments. This 

ability is apparent, except for the variation of E for the previously 

stated reasons, in the error range found for E and G , which is 

x xy 

similar to the standard load cases. 

The overall response was quite good for this target material and 
demonstrated the ability to selectively tailor properties by means of 
the strain energy optimization method. There was no effort made to 
constrain all six parameters. It is felt that this may not be a 
realistic design situation. The failure of the method to return to 
the [0,±45,0] s laminate used for the target properties may be explained 
by the fact that the penality function, with constraint error weights 
included, may be smaller than the strain energy alone for the target 
laminate. The design found by this method may be termed a more 
feasible one if bounds can be formulated for the desired properties. 

5. 1.1. 2 COMPOSITE LAMINATE TAILORED TO MATCH THE PROPERTIES 
OF TITANIUM 

This section deals with tailoring properties to match those of 
the material which is most closely in competition with composite 



TABLE 5 

COMPOSITE LAMINATE - TITANIUM RESULTS FOR PLY ANGLES PERTURBED 
WITH MATERIAL PROPERTIES CONSTRAINED 



Target 

Material 

1 

E ,E and G 
x* y xy 

Constrained 

Cases 

2 

E x and E y 
Constrained 

3 

a x and a y 
Constrained 

Laminate 


[±28,±62] s 

[±77,±13] s 

[±59,±31 ] g 

E (msi) 

* % 

15.80 

9.12 

-42.3 

14.83 

-6.1 

7.66 

-49.2 

E (msi) 
y % 

15.80 

9.12 

-42,3 

14.83 

-6.1 

7.66 

-49.2 

%v( msf ) 

xy % 

6.00 

5.74 

-4.4 

2.28 

-62.1 

6.38 

+6.3 

V 

xy 

% 

0.34 

0.46 

+36.2 

0.13 

-62.7 

0.55 

+61.8 

a (yin/in/°F) 
x % 

4.80 

3.52 

-26.6 

3.52 

-26.6 

3.52 

-26.6 

a (yin/in/°F) 
y % 

4.80 

3.52 

-26.6 

3.52 

-26.6 

3.52 

-26.6 

Stain energy 


1103 

1741 

1061 


NOTE; All percentages refer to variations with respect to the target 
material. 
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laminates from the standpoint of strength, stiffness and weight. Three 
cases were run for this target material, all with standard loadings and 
the following constraints; 

(1) E , E and 6 constrained; 

x y xy 

(2) E and E constrained; 

x y 

(3) a and a, constrained. 

x y 

The results of these examples may be found in Table 5. 

These cases were constrained to have properties equal in both x 

and y directions. This may not be a reasonable situation. The results 

of case (1) show good agreement only for G . The agreement of E , E 

xy x y 

and G simultaneously may not be possible due to the relatively high 
xy 

stiffness in all three directions of titanium. This may be shown in 

case (2) where the constraint on G is relaxed and E and E are more 

xy x y 

readily satisfied. Continuing to case (3), the thermal values are not 

satisfied, but the values determined are equivalent to those of a cross 

ply laminate, thus this is as close as it is physically possible to 

come to the constrained situation. 

This material illustrates the need for selective tailoring. A 

cross-ply laminate [0,90,90,0] s would be a viable solution but only 

by paying a large penalty in strain energy and still violating both 

G and v very seriously. Case (2) best satisfies E and E but 
xy xy x y 

also has the highest strain energy. It may also be noted that each 
of these cases resulted in a laminate which contains two sets of plies, 
each 90 degrees opposed. 

Although results were not as appealing as in the previous section. 
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for a known feasible design, this section reveals the possibility of 
mutually inconsistent design parameters. This highlights the value of 
selective tailoring as an essential in laminate design* 

5. 1.1. 3 COMPOSITE LAMINATE TAILORED TO MATCH THE PROPERTIES 
OF ALUMINUM 

The cases run for aluminum present a different type of problem. 

The elastic and shear modulus for aluminum is low enough that matching 
these properties is not difficult but tne thermal properties are out of 
the feasible range of composite laminates. For this reason the case 
for matching thermal properties was dropped and only the following two 
cases were run for the standard loading conditions; 

(1) E , E and G constrained; 

x y xy 

(2) E and E constrained. 

x y 

As expected these results (Table 6) were good for the stiffness 
properties anl inadequate for the thermal properties. The relaxation 
of the shear constraint resulted in some slight improvement of stiff- 
ness properties but not appreciably. Once again the resultant laminate 
has approximately two sets of plies 90 degrees opposed, but at a 
slightly different orientation. 

This case shows the readily available opportunity for replacement 
of aluminum hardware as long as thermal expansion is not critical. 

5.2 DESIGNS WITH MATERIAL PROPERTIES AND PLY FAILURES CONSTRAINED 

Strength requirements acquire greater importance in primary 



TABLE 6 


COMPOSITE LAMINATE - ALUMINUM RESULTS FOR PLY ANGLES PERTURBED 
WITH MATERIAL PROPERTIES CONSTRAINED 


Cases 

1 2 



Target 

E x’ E y and G xy 

E x and E y 


Material 

Constrained 

Constrained 

Laminate 


[±69»±21] s 

[±64,±28] s 

E (msi) 

10.40 

12.25 

9.29 

x % 


+17.8 

-10.7 

E (msi) 

10.40 

12.25 

9.92 

y % 


+17.8 

-4.6 

G xv (msi) 

3.90 

4.06 

5.50 

xy % 


+4.1 

+41.0 

V vv 

0.33 

0.28 

0.42 

xy 

% 


-15.6 

+27.3 

a (yin/in/°F) 

13.10 

3.52 

3.63 

x % 


-73.1 

-72.3 

a (yin/in/°F) 

13.10 

3.52 

3.42 

y % 


-73.1 

-73.9 

Strain energy 


1276 

1121 


Note: AH percentages refer to variations with respect to the target 
material . 
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structure applications. In these areas, the ultimate load carrying 
capacity is most generally the governing factor in design. The consid- 
eration of a strength criteria in the design is included by means of 
one constraint for each ply failure computation. The Tsai -Hi 1 1 failure 
criteria is used where the numerical value must be less than or equal 
to unity to assure that the individual ply has not failed. The laminate 
is considered failed if one ply has failed. 

On this basis two example problems were run, first by perturbing 
the ply angles for a set number of plies and secondly by perturbing 
the number of plies for a set ply angle with constraints only on ply 
failure. The results of these problems may be found in Table 7. 

These results are for the standard load conditions and the design 
material used previously. The mechanical loads used should evolve an 
initial design such that the worst case is present from a strength 
standpoint. As can be observed from the results, ply failures occur 
for both methods and furthermore, these failures take place in plies 
which are roughly perpendicular to the shear force. The placement of 
these plies is usually the result of the balanced plus-minus angle 
requirements. 

The relative emphasis to be placed on each aspect of the design, 
either property matching or ply failure, must be decided by the 
designer. This emphasis may be implemented by means of load factoring, 
constraint error weighting or some equitable combination of both. The 
examples presented use the same data throughout for the sake of 
comparison. 
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TABLE 7 


COMPOSITE LAMINATE RESULTS FOR PLY ANGLES OR NUMBER OF 
PLIES PERTURBED WITH PLY FAILURES CONSTRAINED 


Cases 

1 2 



Ply 

Number 

Ply Angles 
Perturbed 

Number of Plies 
Perturbed 

Laminate 


[±30,90 ,±21 ] s 

[90, 90, ±45, ±45,0] 


1 

0.6015 

0.6240 


2 

2.1987 

0.6240 


3 

0.7666 

0.2974 

ply failure > 1 

4 

2.0196 

1.0911 


5 

0.7034 

0.2974 


6 


1.0911 


7 

- 

0.4832 

strain energy 


1084 

748 
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5.2.1 DESIGNED BY PERTURBING PLY ANGLES FOR SET NUMBER OF PLIES 

The method of perturbing ply angles for a set number of plies 

relies on the computation of an initial design based on the applied 
loads. This design is then subjected to a series of angular reorienta 
tion until an optimum design is found. From this standpoint the ini- 
tial design provides the number of plies to be considered, within the 
restriction that any plus angle is balanced with a minus angle. 

The loads in the cases run will result in a [0,+45,90] s initial 
laminate design which automatically becomes a [0,±45,90] s laminate due 
to the above restriction. Thus the resultant design could be between 
6 and 12 plies with the factor of two due to symmetry as discussed 
before. 

5.2. 1 .1 COMPOSITE LAMINATE TAILORED TO MATCH THE PROPERTIES OF A 
[0,±45,0] $ LAMINATE 

The all composite material design is geared toward matching a 
known feasible laminate and satisfying ply failure requirements. Com- 
parisons are also made for strain energy in each laminate. Three 
cases were run for the standard load conditions and the following 
constraints; 

(1) E , E and G constrained; 

x y xy 

(2) E and E constrained; 

x y 

(3) a and ct constrained. 

x y 

The results of which can be found in Table 8. Cases (1) and (3) 
yield the best convergence with (3) being the superior design from a 



TABLE 8 


COMPOSITE LAMINATES - [0,±45,0] s RESULTS FOR PLY ANGLES PERTURBED 
WITH MATERIAL PROPERTIES AND PLY FAILURES CONSTRAINED 



Target 

Material 

1 

E ,E and G 
X * y xy 

Constrained 

Cases 

2 

E and E 
x y 

Constrained 

3 

a x and “y 
Constrained 

Laminate 

[0,±45,0] s 

[±29,90,±22] $ 

[+1 ,±89,±15] s 

[±27 ,+51 ,±21] 

E (msi) 
x % 

17.05 

15.43 

-9.5 

19.67 

+15.4 

12.30 

-27.9 

E y ( « si> 

5.23 

8.05 

+53.9 

11.80 

+125.6 

4.76 

-9.0 

V rosi) 

4.43 

4.29 

-3.2 

1.52 

-65.7 

5.73 

+29.3 

v xy 

% 

0.67 

0.44 

-34.3 

0.10 

-85.1 

0.83 

+23.9 

a (pin/ in /°F) 

x % 

2.28 

2.81 

+23.2 

3.11 

+36.4 

2.12 

-7.0 

a (pin/in/°F) 
y % 

6.45 

4.69 

-27.3 

4.18 

-35.2 

6.16 

-4.5 

Strain energy 

% 

1502 

1085 

-27.8 

1590 

+5.8 

923 

-38.6 

Number of plies 


4 

4 

6 


f ai 1 ed 


Note: All percentages refer to variations with respect to the target 
■/ material . 
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property standpoint. Case (2) seems to be invalid and raises the 
possibility that a local minima was found due to the glaring disparity 
of values. Case (1) shows good stiffness tailoring, although the E y 
value is high, which may or may not be beneficial. The thermal 
tailoring characteristics of case (3) are quite good and as a by- 
product give reasonable stiffness values. 

Once again ply failure occurs in the plies perpendicular to the 
applied shear force. This is due to the high value of transverse 
stress in each case. The method is reasonable from the property 
matching standpoint and exhibits capability in both stiffness and 
thermal tailoring in addition to greatly reduced strain energy. 

5.2. 1.2 COMPOSITE LAMINATE TAILORED TO MATCH THE PROPERTIES 
OF TITANIUM 

Three cases were run for standard loads and the constraints that 
follow; 

(1) E v » E w and G vw constrained; 

x y xy > 

(2) E and E constrained; 

x y 

(3) a and a constrained. 

x y 

The results closely parallel those of the preceding section 
(Table 9). The thermal case yielded the best results in all respects 
but all suffer from the same fault discussed in section 5. 1.1. 2 
dealing With the possibility of mutual exclusion of constraint values. 

The method suffers similar ply failure characteristics dealing 
with high transverse tensile stresses in plies perpendicular to the 



TABLE 9 


COMPOSITE LAMINATE - TITANIUM RESULTS FOR PLY ANGLES PERTURBED 
WITH MATERIAL PROPERTIES AND PLY FAILURES CONSTRAINED 




1 

Cases 

2 

3 


Target 

E x> E y a " d G xy 

E x and E y 

a x ancl a y 


Material 

Constrained 

Constrained 

Constrained 

Laminate 


[±55,±74,±43] s 

[±34,±75,±64] s 

[0 j±89,±29] 

E (msi ) 

15.80 

4.36 

6.11 

14.10 

x % 


-72.4 

-61.3 

-10.8 

E (msi ) 

15.80 

11.49 

14.54 

13.77 

y % 


-27.3 

-8.0 

-12.8 

G y „(msi) 

6.00 

5.96 

4.96 

2.95 

xy % 


-0.7 

-17.3 

-50.8 

V vw 

0.34 

0.34 

0.26 

0.18 

xy 

% 


0.0 

-23.5 

-47.1 

a (uin/in/°F) 

4.80 

6.43 

5.46 

3.50 

x % 


+34.0 

+13.8 

-27.1 

a (pin/in/°F) 

4.80 

1.96 

2.48 

3.55 

y % 


-59.2 

-48.3 

-26.0 

Strain energy 


926 

929 

1230 

Number of plies 
failed 


6 

6 

2 


Note: All percentages refer to variations with respect to the target 

material. 
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shear force. The tailoring ability remains quite good, relative to 
the constraint values required. 

5. 2. 1.3 COMPOSITE LAMINATE TAILORED TO MATCH THE PROPERTIES 
OF ALUMINUM 

The cases are, once again restricted to standard loads and con- 
straints on; 

(1) E , E and G constrained; 

x y xy 

(2) E and E constrained. 

x y 

The results yield good stiffness tailoring for each case, and quite 
naturally, poor thermal tailoring (Table 10). Additionally, each case 
yields a laminate not far removed from the initial design. This 
coupled with the ply failure problem would indicate that the initial 
design should contain a larger number of plies. 

5.2.2 DESIGNED BY PERTURBING THE NUMBER OF PLIES FOR A SET PLY ANGLE 

The method presented here relies on the addition or substraction 
of plies for each set ply angle. Ply angles are in five groups con- 
sisting of 90°, ±60°, ±45°, +30° and 0°. Initial designs consist of 
90°, ±45° and 0° angle combinations depending on the applied loads. 

The optimization proceeds by perturbing the number of plies, 
for each set angle, between 0 and 2 until an optimum design is found. 

The strain energy is minimized with respect to the constraints on 
strength and properties. Varying the number of plies should satisfy 



TABLE 10 


COMPOSITE LAMINATE - ALUMINUM RESULTS FOR PLY ANGLES PERTURBED 
WITH MATERIAL PROPERTIES AND PLY FAILURES CONSTRAINED 


Cases 

1 2 

Target E^, E y and G xy E x and E y 



Material 

Constrained 

Constrained 

Laminate 


[0,90,±51] $ 

[±57,90,±5] s 

E (msi ) 

A 0/ 

1° 

10.40 

10.59 

+1.8 

13.77 

+32.4 

E (msi) 
y % 

10.40 

13.18 

+26.7 

12.95 

+24.5 

G xv (msi) 

xy o /q 

3.90 

4.28 

+9.7 

3.35 

-14.1 

v xy 

% 

0.33 

0.27 

-18.2 

0.22 

-33.3 

a (yin/in/°F) 
x % 

13.10 

3.84 

-70.7 

3.45 

-73.7 

a (yin/in/°F) 

y % , 

13.10 

3.26 

-75.1 

3.60 

-72.5 

Strain energy 


1254 

1151 

Number of plies 
failed 


6 

4 


Note: All percentages refer to variations with respect to the target 

material . 
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the strength requirements and provide enough variety of angular orienta- 
tion to satisfy material property constraints. 

5.2. 2.1 COMPOSITE LAMINATE TAILORED TO MATCH THE PROPERTIES OF 
A [0,±45,0] $ LAMINATE 

The cases for this target material use standard load conditions 
with the following constraints; 

(1) E , E and 6 constrained; 

x y xy 

(2) E and E constrained; 

x y 

(3) a and a constrained. 

x y 

These results are summarized in Table 11. Of the three cases only 
(3) responded with viable results. Variations of E have been dis- 
cussed previously and those arguments may also apply here. The 
failure of the stiffness constraints to yield material properties 
of an appropriate nature may be due to the large number of ply 
failure constraints overwhelming the other constraints. 

Ply failure due to transverse stresses of negative angle plies 
continues to appear. The explanation for this may be involved in 
the use of strain energy as the function to be minimized. Another 
possibility that offers an explanation is that of stiffness con- 
straints and ply failure constraints, for this loading case, may be 
mutually exclusive. 

The thermal case offers reasonable tailoring along with a large 
reduction in strain energy. The ply failure constraints for the 
thermal case are close to being satisfied. The total outlook for 



TABLE 11 


COMPOSITE LAMINATE - [0,±45,0] s RESULTS WITH NUMBER OF PLIES 
PERTURBED WITH MATERIAL PROPERTIES AND PLY FAILURES CONSTRAINED 





Cases 




1 

2 

3 


Target 

E x , E y and G 

xy E x and E y 

a „ and a 
x y 


Material 

Constrained 

Constrained 

Constrained 

Laminate 

[0,±45,0] s 

[90,±45,0] s 

[90,90,±60,±45,0] s 

[90, ±30, ±30,0] 

E (msi ) 

17.05 

11.64 

8.02 

15.51 

x % 


-31.7 

-53.0 

-9.0 

E. (msi ) 

5.23 

11.64 

15.00 

7.41 

y % 


+122.6 

+186.8 

+41.7 

G v (msi ) 

4.43 

4.43 

4.43 

4.43 

xy % 


0.0 

0.0 

0.0 

v vu (msi) 

0.67 

0.31 

0.24 

0.49 

xy % 


-53.7 

-64.2 

-26.9 

a (yin/in/°F) 

2.28 

3.52 

4.65 

2.71 

x % 


+54.4 

+103.9 

+18.9 

a (yin/in/°F) 

6.45 

3.52 

2.82 

4.93 

y % 


-45.4 

-56.3 

-23.6 

Stain energy 

1502 

1227 

809 

935 

% 


-18.3 

-46.1 

-37.7 

Number of plies 

6 

4 

6 


failed 


Note: All percentages refer to variations with respect to the target 

material . 
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case (3) Is reasonably good but rather poor for the other two. 

5. 2. 2.2 COMPOSITE LAMINATE TAILORED TO MATCH THE PROPERTIES OF 
TITANIUM 

Three cases were run for titanium with standard loads and the 
following constraints; 

(1) E « E w and G vl/ constrained; 

x y xy 

(2) E and E constrained; 

x y 

(3) ot x and a . 

Results of these cases are found in Table 12 and compare favorably 

with those of the previous section. Cases (2) and (3) have reasonable 

property values but case (1) approaches failure compliance more readily. 

Overall, case (2) would possibly be viable if the E = E constraint 

x y 

were relaxed to some extent, additionally the strain energy is quite low. 

5. 2. 2. 3 COMPOSITE LAMINATE TAILORING TO MATCH THE PROPERTIES OF 
ALUMINUM 

The two cases for aluminum included standard loadings with con- 
straints on the following; 

(1) E y3 E and G constrained; 

x y xy 

(2) E and E constrained. 

x y 

Results found in Table 13 follow the same trend, although the 
stiffness properties are easily match, ply failure of the same nature 
is encountered. 

The aluminum results show that case (1) provides better property 



TABLE 12 


COMPOSITE LAMINATE - TITANIUM RESULTS FOR NUMBER OF PLIES 
PERTURBED WITH MATERIAL PROPERTIES AND PLY FAILURES CONSTRAINED 




1 

Cases 

2 

3 


Target 

E , E and G 
x* y xy 

E x and E y 

a x and a y 


Material 

Constrained 

Constrained 

Constrained 

Laminate 


[90,±60,±60,±30] $ 

[90,±45,±45,0] $ 

[90, 90, ±45,0] 

E (msi) 

A of 

15.80 

6.38 

9.43 

10.18 

(o 


-59.6 

-40.3 

-35.6 

E (msi) 

15.80 

12.63 

9.43 

15.32 

y % 


-20.1 

-40.3 

-3.0 

G xy (msi) 

6.00 

5.43 

-9.5 

5.59 

-6.8 

3.73 

-37.8 

V 

vv# 

0.34 

0.32 

0.45 

0.21 

xy 

% 


-5.9 

+32.4 

-38.2 

a (yin/in/°F) 

4.80 

5.02 

3.52 

4.12 

x % 


+4.6 

-26.7 

-14.2 

a (uin/in/°F) 

4.80 

2.59 

3.52 

3.09 

y % 


-46.0 

-26.7 

-35.6 

Strain energy 


765 

793 

1115 

Number of plies 
failed 


4 

4 

6 


Note: All percentages refer to variations with respect to the target 

material. 



TABLE 13 


COMPOSITE 

LAMINATE - 

ALUMINUM RESULTS FOR 

NUMBER OF PLIES 

PERTURBED WITH MATERIAL 

PROPERTIES AND PLY 

FAILURES CONSTRAINED 




Cases 



1 

2 


Target 

E x> E y and G xy 

E x and E y 


Material 

Constrained 

Constrained 

Laminate 


[90,±45,03 s 

[90,90,±60,±45,0] s 

E (msi) 

10.40 

11.64 

8.02 

x % 


+11.9 

-22.9 

E (msi) 

10.40 

11.64 

15.00 

y % 


+11.9 

+44.2 

G v „(msi) 

3.90 

4.43 

4.43 

xy 2 


+13.6 

+13.6 

v- 

0.33 

0.31 

0.24 

xy 

% 


-6.1 

-27.3 

a {uin/in/°F) 

13.10 

3.52 

4.65 

x % 


-73.1 

-64.5 

a ( pi n/in/°F) 

13.10 

3.52 

2.82 

y 


-73.1 

-78.5 

Strain energy 


1227 

809 


Number of plies 6 4 

failed 


Note: All percentages refer to variations with respect to the target 

material* 
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matching, which is unusual due to extra constraint for shear modulus 
being relaxed In case (2). This may be an example of the strain energy 
function overwhelming the weighted constraint errors and resulting in a 
somewhat misleading design. 



Chapter 6 

SUMMARY AND CONCLUSIONS 

The results and discussion presented in the preceding chapter 
offer several points for observation. 

The results for the optimization without a strength criteria 
proved to be quite good. The ability to tailor a composite laminate 
for desired stiffness or thermal properties was readily apparent. 
Various load conditions provided results which exhibited the respon- 
siveness of the optimization method. Laminates containing a minimum 
of strain energy were found to be feasible to varying degrees, 
dependent upon the constraints applied. 

The inclusion of a strength criteria greatly degraded the quality 
of results. Although the property tailoring capability remained 
viable, it was somewhat restricted. Ply failure, in the plies perpen- 
dicular to the applied shear force was a problem for both angles 
variable and number of plies variable, primarily due to high transverse 
tension stresses. The methods used in this section produced laminates 
with lower strain energy values, also dependent upon constraints 
applied. 

Commenting on the overall design technique, leads to several clear 
points that should come into play during design implementation. The 
first of these deals with selective tailoring of laminate properties. 
Tailoring should be applied on a priority basis with possibly upper 
and lower bounds set for each property. Using a "broad brush" 
approach to tailoring may result in a non-feasible solution situation. 
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Further, along this same vein, constraint values in various directions 
may need to be restricted along with recognizing properties which are 
out of the feasible range of the design material, such as the thermal 
properties of aluminum. Another point that should be noted is that 
of constraint contributions to the penalty function. The effect of 
each constraint should be approximately equal with respect to the 
penalty function. The example problems were all run for the same 
optimization parameters, for comparison purposes. Variation of pro- 
blems would dictate the tailoring of the optimization parameters to 
suit the particular problem. The random point search was found to 
yield excellent starting points for subsequent searches of which the 
creeping search proved the most effective for the example problems. 
Initial designs proved to be somewhat inadequate and should be given 
serious consideration, although this is dependent on applied loads. 

There are some limitations in the design process. The AESOP pro- 
gram is formulated on the basis of an exterior penalty function 
approach to the optimization problem. This means that constraints 
are allowed to oscillate about the desired constraint value instead 
of being bounded by the desired value. This method poses problems, 
especially for failure constraints, where it is essential that the 
values not exceed the desired constraint value. Another limitation 
exists in the number of plies variable method. This method involves 
the use of the number of plies as a continuous design variable. The 
approach is fine for the optimization algorithm but when the function 
evaluations take place the number of plies must be discretized, and 
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thus used in Integer form. Truncation to an integer results in some- 
what erratic behavior in the optimization scheme as was shown in the 
numerical results. 

The first limitation may possibly be alleviated by the use of 
an interior penalty function formulation to impose inequality type 
constraints as opposed to the present equality constraint used. The 
second limitation may be solved by formulating this type of method in 
an integer programming framework, the discussion of which is beyond 
the scope of this investigation. 

The method presented has, in varying degrees, shown that the design 
of a laminate can be accomplished using strain energy minimization as 
the primary criteria. It is felt that this criteria offers a means of 
maximizing the reserve strength that the laminate has the potential to 
exhibit. The inclusion of a failure criteria may be counter to the 
strain energy criteria. This may be viewed by the fact that the energy 
criteria seeks to preserve reserve strength while the failure criteria 
tends to maximize stresses, within failure bounds. Viable trade offs 
may be reached to make these two purposes compatible. 

An optimization method allows the designer the freedom to choose 
directional properties and to eliminate excess material capability in 
non-critical directions. Although, the number of plies variable 
method encountered truncation difficulty, this type of design, using 
families of plies, is the most practical from the present technology 
manufacturing standpoint and as such deserves continued development. 

This type of method coupled with nonlinear programming techniques and 
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continued development of high speed digital computational hardware 
promise to yield a design tool of an extremely powerful nature. The 
use of a tool, such as the AESOP program, should be supplemented with 
both a knowledge of the analysis technique and a background knowledge 
of the optimization procedure. These tools and this knowledge must, 
in the final analysis, be tempered with sound engineering judgment 
which only the experienced designer possesses. 
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A.l LAMAES USER'S GUIDE 

LAMAES is the name chosen for the total design-analysis program, 
consisting of the laminate analysis program coupled with the optimiza- 
tion program, AESOP. 

The LAMAES program is set up to be used in two basic modes, either 
analysis or design. This description is divided into three parts; 
first, information general to both modes, second specifics of the analy- 
sis mode, and third specifics of the design mode. 

A. 1.1 GENERAL 

General data for any computer run, 

CARD 1 format 15 

NUMCAS - number of problems to be run, maximum of 10 
CARD 2 format 1015 

IOPSHN (10) - case mode selection 

0 analysis 

1 design 

Sets of problem data 
CARD 3 format 20A4 

TITLE (20) - title of case 
CARD 4 format 1115 

1. LAMATE - Computation level selection 

0 laminate 

1 lamina 
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2. ISIQEP - lamina input selection 

0 stresses given, strains computed 

1 strains given, stresses computed 

3. I5TRFR - laminate input selection 

1 strains and curvatures given, 
forces and moments computed 

2 forces and moments given, 
strains and curvatures computed 

4. ISYMM - symmetry computation selection 

0 symmetric laminate 

1 asymmetric laminate 

5. IECON - engineering constants computation selection 

LAMATE = 0 LAMATE * 1 

0 none none 

1 laminate only - 

2 lamina and laminate lamina only 

6. NREPT - not used at this time 

7. LAMX - force, moment, strain & curvature computation 
selection 

0 none 

1 all 

2 R - curing thermal loading effects 

3 E - operating thermal loading effects 

note: mechanical effect computed automatically 



LAMY - lamina stresses and strains x,y coordinate 
system computation selection 

0 none 

1 all 

2 M - mechanical loading effects 

3 R 

4 E 

5 M & R 

6 M & E 

7 R & E 

LAMZ - lamina stresses and strains 1»2 coordinate 
system computation selection 

0 none 

1 all 

2 M 

3 R 

4 E 

5 M & R 

6 M & E 

7 R & E 

MAXLAY - option allowing doubling of laminate size 
-1 number of plies remains unchanged 
0 doubles the number of plies symmetrically 
NQPTYP - type of design method option 

-1 no strength criteria considered, angles 
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variable, number of plies constant 

0 strength criteria considered, angles 
variable, number of plies constant 

1 strength criteria considered, number 
of plies variable, for constant ply 
angle 

note: an initial design must be supplied for 

NOPTYP - -1, other options result in 
the automatic determination of an initial 
design (may be unidirectional) 

CARD 5 format 715 

1. IOUTA - material property and layup print option 

1 no printing 

0 printing 

2. IOUTB - lamina stiffness print option 

1 no printing 

0 printing 

3. IOUTC - laminate stiffness print option 

1 no printing 

0 printing 

4. IOUTD - engineering constants print option 

1 no printing 

2 printing 

5. IOUTX - force, moment, strain & curvature* 

6. IOUTY - lamina stress and strain in x,y coordinate 
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system print option* 

7. IOUT Z - lamina stress and strain in 1,2 coordinate 
system print option 

*note : IOUTX, IOUTY and IOUTZ have the 

same combinations, abbreviations 
were previously defined. 

0 none 

1 all 

2 M 

3 R 

4 E 

5 G - total or gross effect 

6 M & R 

7 M & E 

8 M & G 

9 R & E 

10 R & G 

11 E & G 

12 M, R & E 

13 M, R & G 

14 M, E & G 

15 R, E & G 

CARD 6 format 215, 2 F10.0 

NUMLAY - number of plies, maximum 30 

NUMMAT - number of different materials, maximum 30 
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TEMCHG - temperature change, ambient to operating 
DELT 3 - temperature change, cure to ambient 
CARD 7, 8, ... format 7P10.3/5D10.3 

ELP (7, NUMMAT) - lamina material properties 

ELP (1, - ) - elastic modulus E 1 

ELP (2, - ) - elastic modulus E 2 

ELP (3, - ) - shear modulus 

ELP (4, - ) - coefficient thermal expansion, 

ELP (5, - ) - coefficient thermal expansion, a 2 

ELP (6, - ) - Poisson's ratio v^ 2 

ELP (7, - ) - ply thickness 

ELPP (5, NUMMAT) - lamina ultimate strength parameters 

ELPP (1 , - ) - tension ultimate Xy 

ELPP (2, - ) - compression ultimate 

ELPP (3, - ) - tension ultimate Yy 

ELPP (4, - ) - compression ultimate Y^ 

ELPP (5, - ) - shear ultimate Z 
note: there should be 2 cards for each different 

material, 1st card material property and 
2nd card strength parameters. 

CARD 9, ... format 8 (I5,F5.0) see note 
MATER (NUMLAY) - material type 
THETA (NUMLAY) orientation angle 
note: 4 plys per card 



74 


CARD 10 format 6D12.5 
Loads given 

EPSIL (3, 1) for ISIGEP = 1 
or 

SIGMA (3, 1) for ISIGEP = 0 
or 

AMDST (3) 

and for ISTRFR = 1 

AMDCR (3) 
or 

FORCE (3) 
and for ISTRFR = 2 

AMOMN (3) 

A. 1.2. ANALYSIS 

This section deals with specifics of the analysis mode which may 
require various combinations of data values. 

CARD 4 

Use ISIGEP only if LAMATE = 1 
Use ISTRFR only if LAMATE * 0 
Use ISYMM only if LAMATE = 0 

If thermal effects are desired in LAMY then LAMX must 
be compatible, likewise LAMZ must be compatible 
with LAMY. 

Use MAXLAY only if LAMATE = 0, NOPTYP does not apply. 


> 

► for LAMATE = 0 

> 


* 

■ for LAMATE = 1 

i 
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CARD 5 

All print options must be compatible with computations 
or zero's will result. 

CARD 6 

Must supply for analysis. 

CARD 9 

Must supply for analysis. 

CARD 10 

Must supply for analysis. 

A. 1.3 DESIGN 

This section deals with specifics of the design mode which may 
require various combinations of data values. 

CARD 4 

LAMATE = 0 

ISIGEP does not apply 
ISYMM must be compatible with MAXLAY 
IECON * 1 or 2 
LAMX = 1 
LAMY - 1 
LAMZ = 1 
CARD 6 

NUMLAY - 1 or greater even if design is determined internally 
CARD 9 

Must supply for design (1 or greater) 
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CARD 10 

Must supply for design. 

With each set of case data, the control Input for the 
optimization procedure must be input. This is the 
form of FORTRAN NAMELIST type input data. NAMELIST 
input must be preceded by a card denoting the start 
of the NAMELIST input. 

& IAESOP - denotes NAMELIST input to follow, begins in 
column 2. 

NUMOPT - , - number of optimization searches to be 

employed. 

METHOP (20) = , - sequence of searches to be employed 

1 Sectioning 

2 Pattern 

3 Magnify 

4 Steepest Descent 

5 Creeping 

6 Quadratic 

7 Davidon 

8 Random Point 

9 Random Ray 

MAXJJJ - , maximum number function evaluations 

NALPHA = , - number of control parameters to be employed 

ALPHA (100) - , - nominal values of control parameters 

ALPHI (100) * , - upper control parameter search limits 
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ALPLO (TOO) - , - lower control parameter search limits 
NFUNC = , - number of functions to be considered 

NUMPSI = , - number of constraints, maximum of 20 

NPSI (20) = , - constraint function numbers 

2 elastic modulus E direction 

/\ 

3 elastic modulus E y direction 

4 shear modulus 6 direction 

xy 

5 Poisson's ratio v direction 

xy 

6 coefficient of thermal expansion a 

A 

direction 

7 coefficient of thermal expansion a y 
direction 

8-21 ply failure, set internally 
IPRTAL s , - detailed print option for optimization output 

0 no print output 

1 print output for design variables, and 
constraint and objective functions, and 
evaluation number. 

note: many more control variables may be 

set, to obtain more detailed documen- 
tation on the control variables 
see [38]. 

& END 
& IAES0P 
EOF = .TRUE. 
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& END 

This last card ends one complete set of data for a problem, if 
more problems are run then a complete set of data for each would follow 

It may also be noted that the NAMELIST input for a design mode 
problem will override the same input data given in the general data 
section. 

A flowchart (Fig. A.l) describing the basic program branching 
will follow. 
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FIGURE Al. LAMAES 
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Example Input 



i. NUMBER OF CASES 

I 

TEST PROBLEM TITANIUM **' 


0 

2 0 

1 

Oil 

1 0 

0 1 

l 0 

5 

5 5 


1 1 

0 0 

1 

l 


2 1 

LOO. 


-200. 


.308 

.2707 


.9306 .250-5 

.1310-6 

.18806 

-.36206 


.9106 -.6505 

.19205 

1 0. 

1 0. 




.13306 

.13505 

.10206 



C I AESOP 
NFUNC*5, 

NALPHA*2, 

ALPHA* 100*0* f 
IPRTAL*0, 

NUHOPW, 

HETH0?*S,5.2.6 f 3.l,T.13*0, 

MAXJJJ=300, 

NUMPSI*6, 

ALPHI=lOO*90.t 

ALPLO*lO0*-9O., 

5 I BAR* 1.58,1. 58, 4.8, A. 8* 16*. B t 
PS I WT *20*10000.* 

. FTQL*20*.l» 

TTOL»20*l., 

NAXCRP- 5 . 

111*933053127, 

MAXRP T*100, 

CREPHN* 100*1. • 

OCREEP-lOO*!., 

IP*ANL*0, 

CENO 

CIAESe? 

EOF*. TRUE., 

CENO 


21 .0052 



00 

PO 
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Example Output 



*••#•**•**•*•* ***»»»***»*TEST PROBLEM TlTAMlUl 


DESIGN OPTION 

PLY ANGLES PERTURBED - NO STRENGTH CRITERIA 

MATERIAL PROPERTIES CONSTRAINED 
LONGITUDINAL X MOOULUS 
TRANSVERSE V MOOULUS 
THERMAL EXPANSION COEF X 
THERMAL EXPANSION COEF V 


GENERAL COMPUTATION OPTIONS 
NUMBER OF LAYERS 2 
NUMBER OF MATERIALS 1 

TEMPERATURE CHANGE - AMBIENT TO OPERATING 100. 

TEMPERATURE CHANGE - CURE TO AMBIENT -200. 

ANALYSIS LEVEL - LAMINATE « SYMMETRIC 

FORCES AND MOMENTS GIVEN - MtDPLANE STRAINS AND CURVATURES COMPUTEO 


OUTPUT OPTIONS 

LAMX» l 
LAHY» I 
LAM2« l 
10UTA» 0 
IOUTB* 1 
I0UTC* 1 
IOUTO* 0 



INITIAL REPORT 

ALPHA VALUES 

0,00000 0.00000 

FUNCTION VALUES 

1824b. 3.0000 0.27000 2.8000 13.100 0.67407E 06 

**** GAIN *»••* 

HETHOO * 8 OLD PHI * 0.6740716 36 NEW PHI * 0.148282E OS PERCENT GAIN > 97.80317 JJJ > 2 

ALPHA * -0.28056E 02-0. 62959?: 02 

FUNCTN » 0.22054E 04 0.91197E 00 0.9119TE 00 0.352356 01 0.35235E 01 0.148286 OS 

• •**•»»*«•••••*• *««*«•»•• E N o of SEARCH 8 CYCLE NU. I JJJ « 101 

PHI AT START OF SEARCH « 0.674071E 06 PHI AT END OF SEARCH • 0.148282E 05 PERCENT GAIN « 97.80017 


•••••*•«*••****•***•••••• end of Search 5 cycle no. i jjj « H4 **,»***»* Oo 

cn 

PHI AT START OF SEARCH = 0.298165E 05 PHI AT ENO OF SEARCH * 0.26S291E 05 PERCENT GAIN * 10.01926 

END OF SEARCH 2 CYCLE NO. I JJJ * 115 

PHI AT START OF SEARCH * 0.268291E 05 PHI AT ENO OF SEARCH a 0.268291E 05 PERCENT GAIN * 0.00000 


• ••*»•*»•••«•••*•••••••*• ENO OF SEARCH 4 CYCLE NO. I JJJ « 138 

PHI AT START OF SEARCH - 0.515871E 05 PHI AT END OF SEARCH » 0.5138556 05 PERCENT GAIN = 0.39072 


» • 



ENO OF SEARCH A CYCLE NO. 2 Jjj * ZA5 ♦»***•*♦* 

PHI AT START OF SEARCH: * 0.3093666 OT PHI AT ENO OF SEARCH * 3.3 09066E OT PERCENT CAIN » 0.00000 




OF SEARCH 3 

CYCLE NO. 2 

• ♦ ♦ 
JJJ * 

2A8 

* * • r * 
* * * « 

V « » * * 

PHI AT START OF 

SEARCH = 0.61 T753E OT 

PHI AT END OF 

SEARCH « 0.617753E 

07 

* » * 

PERCENT 

* 4 * * 4 

GAIN * 

» * * * • > 

0.00000 



OF SEARCH l 

CYCLE NO. 2 

jjj « 

292 

• * * * 


PHI AT START OF 

SEARCH * 3. 123513E 08 

PHI AT ENO DF 

SEARCH * 0.123S13E 

* * ♦ ♦ * • * • * 9 

08 

* * * 

PERCENT 

* ♦ ♦ i 

CAIN * 

0.00000 




**** 

GAIN 




HETHOO * 

7 OLO PHI « 0.2A6937E 08 

NEW PHI * 

G.2A6987E 08 

PERCENT gain = 

0.00000 

JJJ * 296 

ALPHA * — 0. 1 7*356 

02-0.7J27AE 02 







FUNCTN « 0.291 0 IE OA 0.I369SE 01 0.13695E 01 0.3523SE 01 0.35235E 01 0.2A699E 08 

**** CAIN •*** 

HETHOO « 7 OLO PHI =* 0 .296987E 08 NEW PHI = 0.2A6987E 08 PERCENT GAIN = 0.00000 JJJ = 297 

ALPHA * -0.17A85E 02-0.73269E 02 

FUNCTN = 0.291016 OA 0.13695E 01 0.1369SE 01 0.3S235E 01 0.35235E 01 0.2A699E 08 

»*»« GAIN 

HETHOO * 7 010 PHI • O. 2 A6987E 03 NEW PHI * 0.2AA937E 08 PERCENT GAIN » 0.00000 JJJ * 296 

ALPHA * -0.17A35E 02-0. 73269E 02 


• * • * * • * * * * * • » * * * * • * * * * * * * END OF SEARCH 3 CYCLE -NO. I JJJ ’ l* l • *****«*» 

PHI AT START OF SEARCH = 0.998328E OS PHI AT ENO OF SEARCH = 0.998J28E 05 PERCENT CAIN * 0.30000 

a**************************** ****** a*************************' 


*♦*•*•••**•**•* * : * * END OF SEARCH 1 CYCLE NO. I J JO * 194 **••••*•• 

PHI AT start of SEARCH = 0.196343E 06 PHI AT ENO OF SEARCH = 0.1962036 06 PERCENT GAIN = 0.07124 

**;•■*.*»*;**«••«-*••'•*•*■»•.»■*.**»'*#*• * • »*•*••••*»*»***** *****••*••* I 

a*************************************** OAVIDON SEARCH FAILED TO GA IN**»******«»******»**»»*» **♦•* *********** 

• ••***** •*• ******** • * * * * * END OF SEARCH 7 CYCLE NO. I JJJ * 208 ********* 

PHI AT START OF SEARCH = 0.389590E 06 PHI AT END OF SEARCH * 0.389590E 06 PERCENT GAIN » 0.03000 


a******************* END OF SEARCH 8 CYCLE NO. 2 JJJ = 203 •»•*•**«• 

PHI AT START OF SEARCH = 0.776424E 06 PHI AT END OF SEARCH = V..389590E 06 PERCENT GAIN = 49.82233 


00 


*•«••••• •**«***• **••••••• eno OF SEARCH 5 CYCLE NO. 2 JJJ * 226 **•••***• 

PHI AT START OF SEARCH = 0.154351E 07 PHI AT END 3F SEARCH * 0.154722E 07 PERCENT GAIN * 0.08298 


**••*• • « • • * * ••• *•••**••*• END OF SEARCH 2 CYCLE NO. 2 JJJ = 227 •*••***•* 

PHI AT START OF SEARCH * 0.154722E 07 PHI AT END OF SEARCH * 0.1547226 07 PERCENT GAIN » 0.00000 



•«•* CAIN •*** 


METHOD * 7 OLD PHI = 0. 246937E 03 MEM PHI - 0.246987c 03 PERCENT GAIN * 0.00030 JJJ * 299 

ALPHA * -0.17485E 02-0.73269E 02 

FUHCTN » 0.291016 04 0.13695E 01 0.13695E 01 0.352J5E 01 0.35235E 01 0.24699E 03 

GAIN *««• 

METHOD » 7 OLO PHI * 0.2469876 08 NEW PHI * 0.246987E 03 PERCENT GAIN * 0.00000 JJj =» 300 

ALPHA = -3.17485E 02-0.73269E 02 

FUNCTN ~ 0.291016 04 0. 13695E 01 0.13695E 01 0.3523SE Ct 0.35235E 01 0.24699E 08 

**»• GAIN 

METHOD » 7 OLD PHI * 0.Z46937E 03 NEW PHI » 0.246987E 03 PERCENT GAIN » 0.00000 JJJ » 301 

ALPHA * -0.17485E 02-0.732696 02 

FUNCTN » 0.29101E 04 0.136956 01 0.13695E 01 0.352356 01 0.352356 01 0.246996 08 

FINAL REPORT JJJ * 301 

ALPHA VALUES 
-17.485 -73.269 

FUNCTION VALUES 

2910.1 1.3695 1.3695 3.5235 3.5235 0.24699E 08 



rest PROBLEM TITANIUM 


DESIGN OPTION 

PLY ANGLES P6RTUR3E0 - NO STRENGTH CRITERIA 

HATER I AL PROPERTIES CONSTRAINED 
LONGITUDINAL X ROOULUS 
TRANSVERSE Y MODULUS 
THERMAL EXPANSION COEF X 
THERMAL EXPANSION COEF Y 


GENERAL COMPUTATION OPTIONS 
NUMBER OF LAYERS 3 
NUMBER OF MATERIALS I 

TEMPERATURE CHANGE - AMBIENT TO OPERATING 100. 

TEMPERATURE CHANGE - CURE TO AMBIENT -200. 

ANALYSIS LEVEL - LAMINATE , SYMMETRIC 
NUMBER OF LAYERS INPUT - OOUBLEO 

FORCES AND MOMENTS GIVEN - HIOPLANE STRAINS AND CURVATURES COMPUTED 

CO 

VO 

OUTPUT OPTIONS 

LAMX» 1 
LAMY= l 
LAM2- 1 
10UTA* 0 
IOUTB* 1 
IOUTC* l 
IQUTO 1 0 



LAYER 

NUMBER 

material 

NUMBER 

LATER 

oreintation 

I 

l 

-IT. 

2 

1 

IT. 

1 

1 

-73. 

A 

i 

71. 

5 

;■ l ■ 

73. 

* 

1 

-73. 

t 

1 : 

17. 

« 

l 

-17. 


MATERIAL PROPERTIES 


MATERIAL 

NUMBER 

LONGITUDINAL 

MODULUS 

transver se 

MODULUS 

SHEAR 

MODULUS 

POISSON RATIO 
LONG - TRAN 

POISSON RATIO 
TRAN - LONG 

THERMAL exp 
COEE - LONG 

THERMAL EXP 
COEf - TRAN 

THICKNESS 

l 

0.30000 08 

0.2700D 07 

0.93300 06 

0.21000 00 

0.I893D-01 

0.25030-05 

0.13100-06 

3. 52000-02 


LONGITUDINAL longitudinal transverse transverse 
RATER TENSION COMPRESSION TENSION COMPRESSION 
NUMBER STRENGTH STRENGTH STRENGTH STRENGTH 


SHEAR 

STRENGTH 


l 0.1880D 36 -3.36230 06 3.91000 06 -0.65300 05 0.19Z00 05 


LAYER 

LAMINA COEFFICIENT 

OF THERMAL EXPANSION 

- X.V.2 DIRECT 

NUMBER 

ALPHA X 

ALPHA Y 

ALPHA XY 

l 

0.3*061010-05 

0.12193900-04 

0.59274450-05 

2 

0.34061010-05 

0.12193900-04 

-0.59274450-05 

3 

0.12193900-04 

3.34061010-05 

0.59274450-05 

6 

0.12193900-04 

0. 34061010-05 

-0.59274450-05 

s 

0.12193900-04 

0.34061010-05 

-0.59274450-05 

* 

0. 121.93900-04 

0.34061010-05 

0.59274450-05 

■ ? 

0.34061010-05 

0.12193900-04 

-0.59274450-05 

8 

0.3406 10 I 0-05 

0.12193900-04 

0.59274450-05 


LAMINATE. COEFFICIENT OF THERMAL EXPANSION - X,Y ,1 DIRECTION 
ALPHA X At PHA V ALPHA XY 

0.352353 30-05 0.35235330-05 0*00000000 00 


COMPUTED LAMINATE ENGINEERING CONSTANTS 

LONGITUDINAL TRANSVERSE SHEAR POISSON XY POISSON YX CROSS COEF CROSS COEF CROSS COEF CROSS COEF 

X MODULUS Y MODULUS XY HOOULUS RATIO RATIO XY - X XY - Y X - X» Y - XY 

0.13690 08 0.1369D 08 3.3117D 07 0.19370' 00 0.19370 00 0.00000 00 0.00000 00 O.OOOOD 00 0.0000f> 00 




TOTAL LAMINATE FORCES 

AND DISPLACEMENTS 




FORCE X 

FORCE Y 

FORCE XY 

MOMENT X 

moment y 

MOMENT XY 


3.1191 0440 04 

0.11910440 04 

0.10200000 04 

0.00003000 00 

0.00030000 oc 

0.00300000 

00 

STRAIN X 

STRAIN Y 

STRAIN XY 

Curvature x 

CURVATURE Y 

CURVATURE 

XY 

0.16857150-02 

0.16857150-02 

0.78671210-02 

.0.00000000 00 

3.00303000 30 

0.00000300 

00 



TOrAL LAMINA STRESSES AMD STRAINS 

LAYER 


ARLM3ER 

sigma x 


SIGMA V 


TAU XY 


l 

0.18312770 

0.2 

0.89235650 

09 

0. 1031976D 

05 

2 

0.11131560 

06 

0.18209030 

05 

0.3871870D 

05 

3 

0.89235650 

09 

0.18312770 

02 

0.10319760 

05 

9 

0.18209030 

35 

0.11131560 

06 

0*38718700 

05 

5 

0.18209030 

05 

0. I 113156D 

06 

0.38718700 

05 

■ 6 

0.89235650 

09 

0.18312770 

02 

0.10319760 

05 

7 

0. 11131560 

06 

0.18209030 

05 

0.38718700 

05 

\ ■ 9 '-- : 

0. 18312770 

02 

0.89235650 

09 

0.10319760 

05 


TOTAL LAMINA STRESSES AND STRAINS 


LAYER 

number 

SIGMA l 


SIGMA 2 


TAU 12 


I 

-0.99911900 

09> 

0.13933070 

05 

0.60655690 

09 

2 

0.12500760 

06 

0.95121010 

09 

0.60655890 

09 

3 

-0.999 11900 

09 

0.13933070 

05 

-0.60655890 

09 

9 

0.12500760 

06 

0.95121010 

09 

-0.60655890 

09 

5 

0.12500760 

06 

0.95121010 

09 

-0.60655890 

09 

6 

-0.99911900 

09 

0.13933070 

05 

-0.60655890 

09 

7 

0.12500760 

06 

0.95 1210 10 

09 

0.60655890 

09 

8 

-0.9991 1930 

09 

0.13933070 

05 

0.60655890 

09 


X.Y.2 DIRECTION 


EPSILON X 

EPSILON Y 

GAMMA xy 

0.23669350-02 

0.91299950-02 

0.90526130-02 

0.23669350-02 

0.91299950-02 

0.66816320-02 

0.91299950-02 

0.23669350-02 

0. 93526100-02 

0.91299950-02 

0.23669350-02 

0.66816320-02 

0.91299950-02 

0.23669350-02 

0 .66316320— 02 

0.91299950-02 

0.23669350-02 

0.93526100-02 

0. 23669353-02 

0.91299950-02 

0.66816320-02 

0.23669350-02 

0.91299950-02 

0.90526100-02 


1.2, ■» DIRECTION 

EPSILON l 
-0.13909990-09 
0.93853390-02 
-0.13909990-04 
0.43853390-02 
0.93853390-02 
-0.13909990-09 
0.43853390-02 
— 0.1399999 0-C 9 


EPSILON 2 
0.65G53-»90-J2 
0.21060960-02 
0.65053390-02 
0.21060960-02 
0.21060960-02 
0.65053390-02 
0.21060960-02 
0.65053390-02 


GAMMA 12 
0.65221390-02 
0.65221390-02 
-0 ..65221390-02 
-0.65221390-02 
-0.65221390-02 
-0.65221390-0? 
0.65221390— /2 
0.65221390-02 


LQ 

ro 



STRAIN ENERGY = O.Z910129E 04 


STRENGTH CRITERIA 

LAYER PLY FAILURE 
NUHaCR OCCURS IF >1 

1 2.4*481 

2 0. 17183 

3 2.44481 

4 0. 77183 

5 0.77183 

6 2.44481 

7 0.77183 

8 2.44481 


CD 

OJ 



APPENDIX TABLE A. 3 LAMAES LISTING 


A list of the program is available on request from either author 
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